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RESEARCH  ON  SILICON,  CARBON,  AND 
SILICON  CARBIDE  HETEROSTRUCTURES 

1 .  PROGRAM  OVERVIEW 

The  program  is  a  three-year  effort  to  deposit  thin  films  of 
group  IV  materials  and  heterostructures  in  a  highly  controlled  and 
characterized  deposition  process  and  to  gain  an  understanding  of  the 
process  and  of  the  deposited  materials.  The  deposition  process,  remote 
plasma  chemical  vapor  deposition,  or  RPCVD,  uses  a  flowing  afterglow  of 
metastable  rare  gas  atoms  (in  this  case)  to  dissociate  injected 
reactants  such  as  methane  or  silane  into  desired  deposition  precursors. 
The  process  understanding  is  being  obtained  by  carrying  out  extensive 
characterization  experiments,  such  as  direct  measurement  of  metastable 
densities  and  deposition  precursors,  and  by  modeling  the  kinetics  of 
the  excitation,  recombination,  and  reaction  of  the  gases.  Also  a 
detailed  study  of  the  fundamental  reactions  at  silicon  surfaces  of 
relevant  deposition  precursors,  such  as  the  methyl  radical,  is  being 
carried  out  to  provide  basic  insight  into  the  deposition  process. 
Besides  the  conventional  techniques  that  will  be  used  to  characterize 
the  deposited  films,  a  special  technique,  cryogenic  cathodoluminescence 
measurements  in  an  ultrahigh  vacuum  apparatus,  is  being  applied  to  this 
study . 

The  first  year  of  the  three-year  program  has  been  spent  mainly 
on  assembling  and  characterizing  the  experimental  apparatus.  The 
different  groups  involved  in  the  program  have  been  able  to  proceed 
beyond  this  to  the  stage  of  making  experimental  measurements  relevant 
to  the  program’s  scientific  objectives.  The  four  main  components  of 
the  program,  and  their  respective  technical  leaders  and  locations  are: 
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a.  Thin  film  deposition,  process  characterization,  and 
conventional  materials  characterization:  W.  D.  Partlow, 
Westinghouse  RAD  Center  (Principal  Investigator) . 

b.  Process  modeling:  Dr.  L.  E.  Kline,  Westinghouse  RAD  Center. 

c.  Fundamental  surface  studies:  Prof.  John  T.  Yates  Jr., 
Department  of  Chemistry,  University  of  Pittsburgh. 

d.  Cryogenic  cathodoluminescence  measurements: 

Prof.  W.  J.  Choyke,  Department  of  Physics,  University  of 
Pittsburgh. 

A  technical  overview  of  the  approaches  and  progress  of  all  of 
these  aspects  of  the  program  will  be  presented  in  the  next  four 
sections.  Details  of  the  experimental  and  theoretical  studies  performed 
in  these  program  components  are  reported  in  the  appendices  to  this 
report.  A  final  section  describes  the  publication  and  interaction 
activities  relevant  to  the  program. 
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2.  CHARACTERIZATION  OF  THE 
DEPOSITION  PROCESS 


The  microwave-excited  flowing  afterglow  of  A r  that  is  required 
for  activation  of  deposition  reactions  has  been  experimentally 
characterized.  A  description  of  these  measurements  is  presented  in 
detail  in  Appendix  I.  The  densities  of  the  metastable  states  of  Ar  in 
the  manifold  were  measured  directly  by  excited  state  absorption  of  Ar 
resonance  radiation  and  were  compared  to  the  theoretical  calculations 
that  will  be  described  in  Section  3.  Besides  the  metastable  species 
that  we  measured  in  the  afterglow,  we  found  significant  populations  of 
energy  states,  such  as  the  having  lifetimes  too  short  to  be 

classified  as  metastable.  Several  such  states  were  identified  both  by 
absorption  and  emission  spectroscopy  and  are  expected  to  have  a  major 
influence  on  deposition  chemistry.  We  concluded  that  they  are  at  least 
in  part  due  to  direct  downstream  excitation  of  Ar  by  charged  species, 
since  we  were  able  to  eliminate  their  downstream  optical  emission  by 
inserting  charged  particle-collecting  grids  into  the  flow  stream. 
Experiments  are  continuing  to  fully  control  and  characterize  the 
composition  of  the  afterglow  and  to  determine  the  influence  of  its 
composition  on  the  properties  of  the  deposited  films. 
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3.  MODELING  OF  THE  DEPOSITION  PROCESS 


We  have  modeled  the  microwave-excited  argon  discharge  that 
excites  the  flow  gas  for  creation  of  metastables  in  the  afterglow  for 
the  initiation  of  the  deposition  reactions.  The  operating  point  was 
calculated  for  stable  argon  discharges  at  the  geometry  and  power 
densities  characteristic  of  our  experiments.  This  was  done  by  balancing 
the  production  and  losses  of  electrons  in  the  plasma  by  solving  the 
Boltzmann  equation.  This  determined  the  field/gas  density  ratio,  E/n, 
at  which  the  discharge  operates,  and  which  determines  the  rate 
coefficients  for  the  most  important  processes.  The  appropriate  rate 
coefficients  from  the  literature  were  used  in  a  rate  equation  model  to 
describe  the  plasma  kinetics.  The  resulting  predictions  for  densities 
of  Ar  metastables  were  compared  to  the  experimental  measurements.  The 
predicted  values  have  been  higher  than  the  measured  values,  and  we  are 
examining  the  model  and  critically  evaluating  the  published  data  we  used 
to  explain  the  disagreement.  One  important  quantity,  the  electron 
quenching  rate  for  Ar  metastables,  varies  by  about  a  factor  of  three 
among  four  authors  who  reported  it  in  the  literature.  We  also  are 
taking  into  consideration  the  influence  of  the  high  ionization  fraction 
of  these  particular  plasmas  on  the  shape  of  the  electron  energy 
distribution.  This  may  have  a  dramatic  influence  on  the  predicted 
metastable  densities. 

The  model  predicts  a  dependence  of  metastable  production  on 
power  which  reflects  the  operation  of  the  plasma  in  several  different 
regions  having  different  important  quenching  behavior.  Measured  power 
density  curves  appear  to  reflect  some  of  this  structure. 

These  calculations  are  described  in  more  detail  in  Appendix  II 
of  this  report. 
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4.  CHEMISTRY  OF  REACTIVE  SPECIES 
AT  SILICON  SURFACES 

Construction  and  testing  of  a  UHV  apparatus  for  study  of  the 
reactions  of  hydrocarbon  species  with  Si  surfaces  has  been  completed. 
This  apparatus  has  the  following  features: 

a.  Sample  mounting,  heating,  and  ion  bombardment  capabilities 
needed  for  preparation  and  characterization  of  atomically 
clean  surfaces, 

b.  Absolutely  calibrated  molecular  doser, 

c.  Beam  shutter  for  adsorption/desorption  kinetics  studies, 

d.  Quadrupole  mass  spectrometer, 

e .  Auger  spectrometer . 

A  quantitative  study  of  the  reaction  of  ethylene  with  Si (100)  surfaces 
has  been  carried  out  on  this  apparatus.  The  study,  described  in 
Appendix  III,  was  successful  in  identifying  a  mobile  precursor  state  for 
the  adsorption  of  ethylene  and  determining  the  energetics  for  reaction 
and  desorption  of  this  molecule.  Also,  the  saturation  coverage  and  the 
decomposition  of  desorbed  molecules  were  determined.  A  surface  model  to 
describe  these  data  is  being  developed. 

An  azomethane  source  has  been  prepared.  It  will  be  used  for  the 
study  of  the  reactions  of  methyl  radicals  with  Si  surfaces. 
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5 .  CATHODOLUMINESCENCE  OF 
GROUP  IV  MATERIALS 


The  cathodoluminescence  facility  will  make  it  possible  to  obtain 
characterization  data  highly  relevant  to  the  electronic  properties  of 
the  materials.  This  facility  is  unique  and  takes  advantage  of 
developments  in  vacuum  technology  that  make  it  possible  and  practical  to 
make  good  cathodoluminescence  measurements  at  low  temperatures.  The 
facility  is  now  in  operation,  and  cathodoluminescence  spectra  at 
temperatures  as  low  as  60K  have  been  obtained  on  natural  diamond  samples 
and  on  CVD  diamond  films  that  were  deposited  in  our  laboratory.  Heat 
shields  and  temperature  measurement  techniques  are  presently  being 
improved  to  permit  measurements  near  4.2K.  Details  of  the  facility  are 
described  in  Appendix  IV. 
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6.  PUBLICATIONS  AND  INTERACTION 
ACTIVITIES 

The  publication  activities  are  just  beginning,  as  the 
experimental  results  become  available  on  different  program  tasks.  The 
experimental  work  on  ethylene  interactions  at  Si  surfaces  is  being 
prepared  for  journal  publication.  A  poster  paper  on  the 
cathodoluminescence  facility  -  "An  Ultra  High  Vacuum-Liquid  Helium 
Cathodoluminescence  Monochromator  for  Large  Band  Gap  (Diamond)  Thin  Film 
Analysis,"  W.  J.  Choyke  -  was  presented  at  the  Third  Annual  SDIO/IST-ONR 
Diamond  Technology  Initiative  Symposium  in  July  1988  in  Arlington, 
Virginia. 

Several  publications  that  are  outgrowths  of  the  program  which 
preceded  the  present  one  —  "Plasma  Deposition  of  Silicon  Carbide  Films" 

-  have  appeared  and  will  be  mentioned  because  they  are  relevant  to  the 
program  now  in  place.  A  comprehensive  article  on  the  carbon  film 
deposition  kinetics  work  that  we  concluded  -  "Electron  and  Chemical 
Kinetics  in  Methane  RF  Glow  Discharge  Deposition  Plasmas,"  L.  E.  Kline 
and  W.  D.  Partlow  —  is  scheduled  to  appear  in  the  Journal  of  Applied 
Physics  on  December  15,  1988.  A  collaboration  with  NASA  workers  on 
Silicon  Carbide  films  that  were  developed  on  the  program  -  "Grazing 
Incidence  Reflectance  of  SiC  Films  Produced  By  Plasma-Assisted  Chemical 
Vapor  Deposition,"  Ritva  A.  M.  Keski-Kuha,  John  F.  Osantowski, 

Albert  R.  Toth,  and  William  D.  Partlow  -  was  published  in  Applied  Optics 
27  No.  8,  1499  (April  1988).  A  paper  describing  the  silicon  carbide 
deposition  kinetics  -  "Electron  and  Chemical  Kinetics  in  Low  Pressure  RF 
Discharge  Deposition  of  Amorphous  Silicon-Carbide  Thin  Films," 

L.  E.  Kline  and  W.  D.  Partlow  -  has  appeared  in  the  Proceedings  of  the 
Eighth  International  Symposium  on  Plasma  Chemistry,  Tokyo,  Japan,  August 
1987. 
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APPENDIX  I 

CHARACTERIZATION  OF  THE  DEPOSITION  PROCESS 
I.  Experimental  Apparatus 

The  deposition  reactor  consists  of  a  microwave-excited  flow  tube 
from  which  excited  argon  atoms  dissociate  reactants,  such  as  methane  for 
carbon-containing  films  or  silane  for  silicon-containing  films,  and  the 
reactive  species  thus  produced  are  flowed  over  a  sample  susceptor.  The 
reactor  is  shown  schematically  in  Figure  1  in  side  view,  which  also 
shows  how  grids  are  placed  in  the  flow  stream  to  control  charged 
species.  The  excitation  source  is  a  0-1.5  kW  microwave  source  at  2.455 
GHz  equipped  with  constant  frequency  power  control  and  four-phase  load 
impedance  characterization.  A  four-stub  tuner  is  used  to  match  the 
microwave  excitation  into  a  2.54  cm  cylindrical  plasma  tube  exciter, 
which  has  8.5  cm  cylindrical  extensions  to  attenuate  the  microwave 
fields.  The  flow  tube  expands  to  4.8  cm  diameter  immediately  downstream 
of  the  plasma  excitation  unit.  The  process  pumping  is  done  with  a 
mechanical  pump-backed  Roots  blower  that  has  a  nominal  pimping  speed  of 
1000  liter/second.  The  assembled  system  with  valves  and  downstream 
pressure  controller  in  place  delivers  411  liter/second  pumping  capacity. 
This  will  produce  a  plug  flow  velocity  of  about  2.3(10)  cm/sec  in  the 
4.8  cm  flow  tube.  This  is  adequate  for  transporting  rare  gas 
metastables  from  the  excitation  region  to  the  reaction  zone  without 
significant  losses.  The  sample  susceptor  accepts  several  different 
heating  and  cooling  modules  and  can  provide  temperatures  from  room 
temperature  to  900*C.  It  is  interchangeable  with  a  differentially 
pumped  mass  spectrometer  for  sampling  neutral  and  ionized  species  from 
the  deposition  gases. 
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Figure  1.  Schematic  of  the  remote  plasma  CVD  reactor. 


The  purpose  of  these  characterization  measurements  is  to  obtain 
the  information  needed  to  gain  an  understanding  of  the  production  of 
reactive  species  that  lead  to  deposition,  and  to  control  these  species 
to  change  the  properties  of  the  deposited  films.  The  first  step  is  to 
determine  which  species  from  the  excitation  region  reach  the  reaction 
zone.  We  are  currently  obtaining  this  information  for  argon  flow  gases 
by  using  two  diagnostic  techniques:  optical  absorption  and  emission 
spectroscopy.  A  top  view  schematic  of  the  deposition  reactor  provided 
in  Figure  2  shows  how  these  diagnostics  are  arranged.  The  emission 
spectrometer  consists  of  a  0.25  meter  Jarrell  Ash  monochromator  with 
computer  control.  A  0.3  meter  McPherson  monochromator  is  used  for 
absorption  of  resonance  emission  lines  from  an  Ar  discharge  tube.  The 
light  from  the  absorption  source  is  chopped  at  200  Hz,  and  synchronous 
detection  provides  discrimination  from  any  light  produced  in  the 
deposition  reactor  from  the  excitation  plasma  or  from  species  that  flow 
to  the  reaction  zone.  A  beam  splitter  and  filtered  detector  are  used  to 
monitor  the  output  of  the  resonance  radiation  used  for  the  absorption 
measurements  so  that  some  of  the  advantages  of  double-beam  spectroscopy 
can  be  obtained. 

A  photograph  is  shown  in  Figure  3  of  the  deposition  facility  and 
its  diagnostic  equipment. 


XX.  Optical  Emission  Measurements 

Argon  metastables  in  the  **P  manifold  located  11.5  to  12  eV  above 
the  ground  state  do  not  lead  to  visible  emission,  so  the  visible 
emission  seen  in  the  afterglow  when  no  grids  are  used  is  an  indication 
that  other  excited  species  are  present.  A  photograph  of  the  afterglow 
under  typical  operating  conditions,  presented  in  Figure  4,  shows  this 
emission  clearly,  and  its  visible  emission  spectrum  is  shown  in 
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Figure  2.  Schematic  top  view  of  the  optical  diagnostics  on  the  RPCVD  reactor. 


Figure  3.  Remote  plasma  CVD  experimental  facility. 
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Figure  4.  Photograph  of  flowing  argon  afterglow  w  h  no  grids  in  the 
flow  tube. 
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Figure  5.  At  lov  excitation  powers,  emission  from  the  afterglow 
disappears  when  grounded  grids  are  inserted  into  the  flow  stream, 
indicating  that  the  emission  is  due  to  charged  species  directly  exciting 
the  gas  in  the  afterglow  region.  Some  visible  afterglow  is  seen  at 
higher  powers,  but  there  are  indications  that  this  is  due  to  leakage 
around  the  grids.  We  have  redesigned  the  grids  to  eliminate  this 
leakage  and  plan  to  measure  grid  currents  as  a  function  of  grid  bias, 
and  to  correlate  these  measurements  with  non-metastable  species  detected 
downstream. 


Fiou. i;  /Vjl  LP* i 
PtHAG#.:  4£o  vaJ 
rio  <5*/DS 


Figure  5.  Spectrum  of  the  argon  flowing  afterglow  with  no  grids  in  flow 
tube. 
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1X1.  Optical  Absorption  Measurements 

3 

Excited  state  absorption  from  the  metastable  P  levels  is  easily 
detected  by  the  absorption  diagnostic.  Absorbances  of  over  60%  are 
measured  at  the  highest  powers  at  811.5  am,  the  strongest  absorbing  line 
from  the  state.  An  example  of  the  raw  absorbance  data  is  shown  in 
Figure  6,  where  it  is  plotted  as  a  function  of  pressure.  The  density  of 
the  metastables  can  be  obtained  from  this  absorption  data  with  a 
considerable  effort  of  data  analysis.  The  difficulty  is  due  to  the  fact 
that  the  (Doppler)  linewidth  of  the  Ar  lines  used  to  measure  the 
absorption  is  broader  than  the  linewidth  of  the  measured  Ar  absorption 
line,  and  these  linewidths  are  not  accurately  characterized.  They  also 
depend  upon  temperature.  The  absorption  does  not  obey  the  Beer-Lambert 
absorption  law,  absorbance  values  saturating  at  values  below  unity.  We 
used  a  technique  developed  by  Piper  et  al.,  which  utilizes  absorption 
measurements  made  with  lines  from  the  same  state  having  different 
oscillator  strengths.  We  also  made  measurements  with  different 
metastable  densities.  This  provides  an  effective  in-situ  technique  for 
determining  the  linewidth  information  that  is  required  for  determining 
densities.  It  basically  consists  of  determining  the  kl  values 
(absorption  coefficient  multiplied  by  absorption  length)  for  several 
values  of  a,  the  ratio  of  source  to  sample  linewidths.  Sample  density 
depends  on: 

kl/f, 

where  f  is  the  oscillator  strength  of  the  transition.  The  same  sample 
density  should  be  predicted  from  measurement  of  different  absorption 
lines  from  the  same  state.  The  value  of  a  is  used,  which  gives  the  most 
consistent  density  predictions  for  the  same  density  and  does  this  over  a 
wide  range  of  densities.  The  variation  in  density  was  obtained  by 
varying  the  discharge  power,  and  the  resulting  kl/f  values  as  a  function 
of  power  for  four  different  a  values  are  shown  in  Figure  7.  An 
examination  of  the  data  shows  that  a  values  of  2.5  give  the  tightest 
patterns  over  the  entire  range  of  discharge  powers.  It  is  significant 
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Figure  6.  Absorption  at  811.5  n>  from  the  A r  state  versus  pressure. 


is 


that  lower  a  values  are  not  best  at  high  powers,  as  would  be  expected  if 
there  were  broadening  of  the  afterglow  lines  due  to  heating  from  the 
discharge.  An  additional  confirmation  of  our  selection  of  a  is  obtained 
from  Figure  8,  where  kl/f  versus  f  is  plotted  for  a  single  value  of 
discharge  power.  Systematic  departures  from  the  ideal  independence  of 
kl/f  on  f  are  seen  for  a  values  above  and  below  2.5. 

3 

The  densities  determined  by  absorption  spectroscopy  of  the  P 
levels  at  2  torr  pressure  (without  grids)  are  plotted  in  Figure  9.  The 

ratio  of  the  ^P0  to  the  densities  are  comparable  to  measurements  of 

1  z  u  3 

other  workers.  '  The  significantly  high  density  of  the  state, 

which  has  a  lifetime  of  about  10  microseconds  and  is  not  metastable, 

gives  confirmation  to  our  earlier  conclusion  that  there  is  direct 

excitation  in  the  afterglow.  The  departure  from  linearity  of  the  power 

dependencies  of  the  metastable  densities  is  outside  the  experimental 

error  that  we  ascribe  to  our  absorption  measurements.  The  modeling  work 

to  be  discussed  in  the  next  section  predicts  such  nonlinear  behavior. 


XV.  Future  Plans 

More  experiments  will  be  done  to  determine  the  nature  of  the 
downstream  excitation  and  to  achieve  the  ability  to  control  it.  As  we 
mentioned,  biased  grid  measurements  are  planned  for  this  purpose.  In 
addition,  we  plan  to  test  our  conclusions  by  measuring  emission  from 
methane  and  nitrogen,  both  of  which  are  sensitive  to  the  composition  of 
the  afterglow,  as  well  as  from  silane.  These  measurements,  in 
conjunction  with  the  modeling,  will  put  us  in  a  position  to  relate  the 
properties  of  films  deposited  from  the  afterglow  to  the  state  of  the 
afterglow  excitation. 
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APPENDIX  II 

MODELING  OF  AEGON  DISCHAEGE 
Determination  of  Operating  Point 

The  first  step  in  modeling  the  argon  microwave  discharge  is  the 
determination  of  its  operating  point.  This  means  finding  the  value  of 
the  ratio  of  the  electric  field  to  the  gas  density,  E/N,  at  which  the 
discharge  operates.  Let  (E/N)q  denote  this  value  of  E/N.  Finding 
(E/N)g  is  of  critical  importance  because  it  determines  the  rate 
coefficients  for  many  important  processes,  such  as  the  creation  of 
metastable  states  of  argon. 

The  discharge  will  operate  at  a  steady  state  such  that  the  production 
and  loss  processes  for  electrons  are  in  balance.  The  only  important 
production  process  is  electron  impact  ionisation  of  ground  state  argon 
atoms.  Two-step  ionisation  —  that  is,  electron  impact  ionisation  of 
argon  metastables  —  should  not  be  important  for  expected  metastable 
concentrations.  At  the  anticipated  operating  pressures  of  a  few  Torr  or 
less,  the  dominant  electron  loss  process  should  be  ambipolar  diffusion 
to  the  walls  of  the  discharge  tube.  Volume  recombination  and  free 
electron  diffusion  are  not  significant.  Such  diffusion-dominated 
discharges  are  discussed,  for  example,  by  Cherrington. *  In  a  steady- 
state,  diffusion-dominated  discharge,  the  electron  conservation  equation 
is 


dn  /dt  =  0  =  k.  n  N  -  D  n  /A2, 

where  ng  is  electron  density,  is  the  rate  coefficient  for  electron 
impact  ionisation,  N  is  the  argon  atom  density,  D&  is  the  ambipolar 


»i 


diffusion  coefficient,  and  A  =  (tube  radius) /2. 405  is  the  characteristic 
ambipolar  diffusion  length.  Thus,  in  steady-state  the  ionization  rate 
per  electron,  iA  ,  is  given  by 


v. 

x 


The  two  quantities  to  be  determined  as  a  function  of  E/N  are  and  D&. 

k.  is  found  in  a  straightforward  way  by  solving  the  Boltzmann  equation, 
i  1 

D  is  given  by: 


D  = 
a 


D 

e 


u  +  u 
re  r+ 


In  discharges  such  as  this  where  the  electrons  are  strongly  heated,  it 
is  true  that  D  »  D  and  p  »  p  .  Thus,  we  have  approximately 

Da  =  GW.  K  *  D+  • 

2 

Using  a  literature  value  for  /i+,  we  have 
D  p  =  (D/p  )  (1200)  cm^  Torr/s 

ft  G 

for  argon,  where  p  is  pressure  in  Torr. 

o 

Plots  of  and  D&/A  2  versus  E/N  are  shown  for  0.3  Torr  and  2.0  Torr 

in  Figures  1  and  2,  respectively.  The  predicted  (E/N)r  for  these  two 

-17  ^  2 

pressures  are  60  Td  and  28  Td,  where  1  Td  =  1  x  10  V  cm  . 


aa 


IONIZ./LOSS  RATE  PER  ELECTRON 


0  1  0  20  30  40  50  60  70  80 


E/N,  Td 


Figure  1.  Ionization  rate  and  diffusion  loss  rate  versus  B/N  for  argon 
in  a  cylindrical  tube  of  1.1  cin  radius  at  a  pressure  of  0.3 
Torr.  The  intersection  of  the  curves  is  the  predicted 
operating  point  of  the  discharge. 
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Figure  2.  Ionization  rate  and  diffusion  loss  rate  versus  E/N  for  argon 
in  a  cylindrical  tube  of  1.1  cm  radius  at  a  pressure  of  2.0 
Torr.  The  intersection  of  the  curves  is  the  predicted 
operating  point  of  the  discharge. 


Solutions  of  the  Boltzmann  Equation 


The  cross  sections  used  for  the  4p  states  of  argon  are  those  of 
3 

Tachibana.  There  are  four  4p  states,  two  of  which  are  metastable.  The 
metastable  states  are  the  and  ^Pn  states  at  11.55  eV  and  11.72  eV. 

i  3U  1 

The  two  resonant  4p  states  are  the  at  11.62  eV  and  the  P^  at  11.83 

eV.  The  other  two  cross  sections  found  by  Tachibana  are  for  lumped 

states:  the  higher  lying  optically  forbidden  (mostly  p)  states  and  the 

higher  lying  optically  allowed  s  and  d  states.  The  ionization  cross 

4 

section  is  that  of  Rapp  and  Englander-Golden. 

The  4p  states,  in  addition  to  population  through  direct  excitation,  are 
also  populated  by  cascading  from  the  higher  forbidden  states.  From 
measurements  and  Boltzmann  calculations,  Tachibana  concludes  that  39%  of 
the  total  excitation  of  the  higher  forbidden  states  cascade  to  Pg, 
while  only  7%  cascades  to  ^Pq.  The  remainder  of  these  states  cascades 
to  the  resonant  ^P^  and  ^P^  states. 

Ve  have  performed  Boltzmann  equation  solutions  for  argon  using  the  cross 

sections  of  Tachibana.  We  have  determined  excitation  rate  coefficients 
3 

in  cm  /s  for  0.3  Torr  and  2.0  Torr  as  follows: 


0.3  Torr  (60  Tdl 

2.0  Torr  (28 

% 

5.40(-ll) 

1 . 19 (—11) 

5.59(-ll) 

1.63(-11) 

1.02(-11) 

2.26(-12) 

lpl 

8.63(-ll) 

2.29(-ll) 

3ft 


Equation  Discharge  Model 


The  argon  discharge  has  been  modeled  using  a  rate  equation  computer  code 

developed  by  Lacina.^  The  processes  included  in  the  model  are  shown  in 

Table  1  along  with  their  reaction  rates.  Processes  1-4  are  electron 

impact  excitation  of  the  two  metastable  and  two  resonant  4p  states.  The 

rates  for  these  are  given  above.  Processes  5  and  6  are  electronic 

quenching  of  the  metastables.  Quenching  rates  have  been  given  by  four 

authors  and  vary  by  about  a  factor  of  three.  The  rates  adopted  here  are 

6  3..  —7  3 

those  of  Bochkova  and  Sukiasyan  for  the  °Pn  state  (7  x  10  cm°/s)  and 

7  3  °  7  3 

Gerasimov  and  Petrov  for  the  state  (4.5  x  10  cm  /s) . 

Reaction  7  is  ionization  with  the  rate  from  our  Boltzmann  calculations. 
Reactions  8  and  9  are  two-  and  three-body  quenching  of  metastables, 
which  are  included  although  they  have  not  had  a  significant  effect  on 
any  of  the  calculations  thus  far.  The  rates  are  taken  from  Kolts  and 

g 

Setser.  Reactions  10  and  11  correspond  to  spontaneous  radiation  from 
the  resonant  states.  The  rates  include  the  effect  of  radiation  trapping 
and  are  estimates  from  Scheller  et  al.  Reaction  12  represents  the  loss 
of  positive  ions  due  to  ambipolar  diffusion  to  the  tube  wall.  It  is 
treated  as  a  volumetric  loss  with  a  rate  determined  by  the  discharge 
operating  point. 

The  electron  density  is  held  fixed  at  a  value  chosen  to  match  the 
expected  discharge  power  density 

E  j  =  (E/N) N  n  e  W  . 

Here  e  is  the  electronic  charge  and  W  is  the  electron  drift  velocity. 

E/N  and  V  are  determined  by  the  discharge  operating  point  through 
pressure  and  tube  radius.  Thus,  according  to  the  model,  only  ng  changes 
to  correspond  to  varying  experimental  discharge  power  density. 
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Note:  Reactions  no.  7  and  no.  12  appear  charge-unbalanced  because  of  a  formatting 


We  simulated  experimental  conditions  of  2.0  Ton*  and  a  power  density  of 

3  11 

about  2  watts/cm  ,  corresponding  to  an  electron  density  of  2.5  x  10A1 

-3 

cm  .  The  resulting  steady-state  densities  are  shown  in  Table  2.  The 
3  12-3 

Pg  density  of  6.45  x  10  cm  is  about  40  times  higher  than  the 
measured  value.  Figure  3  shows  the  3P2  density  as  a  function  of  time 
for  the  first  500  microseconds  of  discharge  operation.  In  gas  flowing 
at  250  m/s,  this  corresponds  to  a  flow  distance  of  12.5  cm  so  the 
metastable  density  attains  its  full  value  in  about  2  cm. 

Table  2.  Predicted  steady-state  species  densities  of  an  argon  discharge 

with  a  pressure  of  2  Torr  and  electron  density  of  2.5  x  1011 

-3  q 

cm  .  The  ar(m)  and  ar (m2)  species  are  the  *r0  and  r_ 

«  u 

metastable  states,  while  the  ar(r)  and  ar(r2)  are  the  resonant 
states  and  ar(p)  is  the  positive  ion.  Excitation  energy,  rate 
of  population  change,  and  time  constant  are  also  given  for 
each  state. 


time  t  -  5.000E-04  sec 


population  densities  and  rates  of  change 


i 

species 

name 

e<i), 

(ev) 

n(i) 
(cm-3 ) 

dn( i )/dt 
(cm-3/sec) 

tau(  1 
1.0E-06  si 

1 

rad 

0.00 

0 .  OOOE+OO 

0. 000E+00 

2 

e(-) 

0.00 

2 . 500E+11 

0. 000E+00 

3 

ar 

0.00 

6 . 681E+16 

2.211E+13 

3.02E+09 

4 

ar(m) 

11.55 

6 . 454E+12 

-8.151E+12 

7.92E+05 

5 

ar(n>2) 

11.62 

4 . 418E+12 

-2. 843E+12 

1.55E+06 

6 

ar(r) 

11.71 

4 . 698E+12 

-5 . 869E+12 

8.00E+05 

7 

ar(r2) 

11.83 

4 . 203E+12 

-2 . 147E+12 

1 . 96E+06 

8 

ar(p) 

15.80 

2 . 369E+11 

-3 . 102E+12 

7.64E+04 

total  gas  density 

6 . 683E+16 

integration  step  size  =  69.269 

order  of  integration  =  4 
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Dependence  of  Metas table 
Production  on  Power  Density 


The  only  significant  process  for  metastable  destruction  is  electron 
quenching,  which  produces  a  resonant  state  that  can  then  radiate  to  the 
ground  state.  The  estimated  effective  lifetime  of  the  resonant  states 
is  10  microseconds  due  to  radiation  trapping.  This  radiative  lifetime 
competes  with  a  lifetime  due  to  electron  impact  of  the  resonant  state 
that  produces  the  metastable  state.  The  rate  for  this  process  is  given 
by 


"ii  -  \  % 


where  k^  is  the  rate  coefficient.  Thus,  for  sufficiently  high  electron 
density,  this  process  can  dominate  radiative  decay  and  electron 
quenching  is  defeated.  But  in  the  other  limit  of  low  electron  density, 
all  quenched  metastables  are  permanently  lost  because  of  radiative 
decay.  In  this  case,  the  ratio  of  metastable  to  ground  state  densities 
is  given  by 


[Ar*]/[Ar]  =  k+  / 


where  k  is  the  rate  coefficient  for  production  (reaction  1)  and  kn 

*  .  .  3_  #  *  % 

is  the  rate  coefficient  for  rg  quenching  (the  forward  rate  for  reaction 
5).  Using  the  rate  coefficients  from  Table  1,  the  fractional  metastable 
density  is  2.64  x  10  For  a  pressure  of  2  Torr  this  gives  a  predicted 
metastable  density  of  1.77  x  10  cm  in  the  limit  of  low  electron 
density. 


For  the  case  shown  in  Table  2  and  Figure  3,  the  rate  for  electron  impact 
de-excitation  (1.13  x  10^  s  *  is  nearly  equal  to  the  rate  for  radiation 

(10  s  ) .  In  this  case  the  r„  density  is  about  four  times  the  low  n 

*  12  ® 
value.  Increasing  the  electron  density  by  a  factor  of  four  to  10A 


39 


cm  ^  increases  the  density  by  a  factor  of  3.2.  A  further  order  of 
magnitude  increase  results  in  the  ^Pg  density  rising  by  a  factor  of  9.2. 
Thus,  the  model  predicts  three  regions: 

1)  For  low  electron  density,  metastable  density  is  independent  of 
electron  density  (power  density) . 

2)  For  "intermediate*  electron  density,  i.e.,  where  the  radiative  and 

3  1 

collisional  depopulation  of  the  resonant  P^  and  P^  states  are 
comparable,  there  is  a  transition  region  where  metastable  density 
increases  sublinearly  with  electron  density. 

3)  For  high  electron  density  the  metastable  density  increases  linearly 
with  electron  density. 
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population  density  n  (cra-3)  and  rate  of  change  dn/dt  (cm-3/s)  for  ar(m) 
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Figure  3.  Predicted  spatial  behavior  of  the  argon  metastable  density  (***)  and  its  tine  derivative 
(+++).  The  horizontal  axis,  which  is  labelled  "time*  and  goes  to  500  microseconds, 
corresponds  to  a  total  distance  of  12.5  cm  at  the  flow  velocity  of  250  m/s. 


APPENDIX  HI 

THE  REACTION  OF  ETHYLENE  WITH  Si  (100) 


I.  INTRODUCTION 

Understanding  the  interaction  of  simple  hydrocarbon  molecu¬ 
les  with  silicon  single  crystal  surfaces  is  of  importance  in 
determining  the  mechanism  of  formation  of  silicon  carbide  films 
on  Si.  The  interaction  of  hydrocarbons  with  the  dangling  bonds 
at  the  Si  surface  represents  the  very  first  step  in  the  silicon 
carbide  thin  film  formation  process.  We  have  employed  surface 
kinetic  and  surface  spectroscopic  methods  to  study  the  details  of 
adsorption,  thermal  decomposition  and  reaction  of  the  simplest 
olefin,  ethylene,  on  the  Si(100)  surface. 

It  has  previously  been  shown  in  our  laboratory  [1]  and 
elsewhere  [2]  that  the  carbon-carbon  double  bond  is  an  active 
center  for  reaction  with  atomically  clean  Si  surfaces. 

Comparison  experiments  involving  alkane  and  olefin  molecules 
clearly  show  that  at  120  K  only  the  olefins  react  efficiently 
with  the  Si  dangling  bonds  on  Si (100)  [3]. 


II.  THE  Si (100)  SURFACE 

It  is  well  known  that  the  Si (100)  surface  consists  of  rows 
of  Si-Si  dimers  exhibiting  a  (2  x  1)  lattice  structure  [4].  Each 
of  the  Si  atoms  in  a  dimer  species  contains  an  upward-directed  Si 
dangling  bond  which  can  act  as  a  reactive  center  for  the  adsorp¬ 
tion  of  olefin  molecules.  A  diagram  of  the  Si  dimer  sites  on 
Si(100)  is  shown  in  figure  1. 
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Figure  1.  Ideal  Si(100)(2xl)  surface.  The  open  and 
hatched  circles  represent  the  atoms  in  the 
topmost  and  second  layers,  respectively.  (a) 
Unreconstructed  surface.  (b)  Reconstructed  sur 
face  containing  Si-Si  dimer  species  [5].  There 
are  3.4  x  10^4  si  dimers/cm^. 


III.  EXPERIMENTAL  METHODS 

A.  Preparation  of  the  Si (100)  Surface 

The  Si (100)  crystal  (10  ohm-cm,  B  doped)  is  held  between  W 
heating  support  wires  using  Ta  sheet  spring  clamps  which  grip  the 
crystal  internally  in  slots  cut  along  the  edges  of  the  rec¬ 
tangular  crystal  [6].  In  addition,  a  W/5%Re :W/26%Re  thermocouple 
junction  is  contained  inside  a  third  slot  in  the  crystal  edge. 
This  arrangement  permits  ohmic  heating  of  the  crystal  for 
cleaning  as  well  as  for  temperature  programmed  desorption 
spectroscopy  measurements  [6,7].  The  crystal  is  ion  bombarded 
with  Ar+  (2k  eV;  2.5y  A/cm^;  1800  seconds)  for  initial  cleaning 
and  is  then  heated  to  1190  K  to  remove  the  last  traces  of  surface 
carbon  and  oxygen  as  measured  by  Auger  spectroscopy.  Care  must 
be  taken  not  to  heat  the  crystal  to  temperatures  above  1350  K  in 
order  to  avoid  irreversible  damage  to  the  crystal  as  observed  by 
STM  measurements  [8].  Following  cleaning,  the  crystal  is  cooled 
to  about  100  K  for  adsorption  studies. 
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B.  Use  of  a  Calibrated  Directed  Beam  Doser  for  Adsorption 

Fundamental  insight  into  the  stoichiometry  of  adsorption  on 
single  crystal  surfaces  may  be  made  by  employing  quantitative 
adsorption  measurements.  This  may  be  achieved  through  the  use  of 
a  directed  beam  doser  of  our  own  design  [6]  which  permits  the 
establishment  of  an  absolute  known  flux  of  adsorbate  molecules 
which  impinge  on  the  single  crystal  surface.  Through  the  use  of 
a  mass  spectrometer  for  the  detection  of  the  reflected  flux  of 
molecules,  it  is  possible  to  determine  both  the  absolute  coverage 
and  the  absolute  reaction  probability  of  the  impinging  molecules 
with  the  surface.  The  angular  distribution  of  the  gas  from 
dosers  of  this  design  has  been  carefully  investigated  [9]. 


C.  Apparatus  for  Adsorption-Desorption  Kinetic  Studies 

Figure  2  shows  a  schematic  diagram  of  the  apparatus  which 
has  just  been  built  for  the  study  of  hydrocarbon  molecule  adsorp¬ 
tion  on  Si(100).  It  consists  of  a  beam  doser  which  delivers  a 
known  flux,  FC2H4'  ethylene  to  the  crystal.  The  apparatus  is 
equipped  with  a  shielded  quadrupole  mass  spectrometer  which  is 
arranged  to  receive  only  a  random  flux  of  gas  which  is  either  not 
adsorbed  by  the  crystal  or  which  misses  the  crystal  and  rebounds 
from  the  walls  of  the  pumped  ultrahigh  vacuum  chamber.  A  shutter 
may  be  rapidly  moved  into  the  C2H4  beam  to  protect  the  crystal 
from  direct  adsorption  from  the  beam.  The  lower  portion  of 
figure  2  shows  a  schematic  experiment  in  which  adsorption  by  the 
Si (100)  crystal  is  monitored  by  the  mass  spectrometer.  A 
background  level  of  C2H4  is  initially  present  in  the  apparatus? 
the  beam  doser  is  turned  on  to  admit  C2H4  with  the  shutter  in 
place  to  intercept  the  direct  beam,  and  after  a  short  time,  a 
stable  pressure  of  C2H4  is  observed  by  the  mass  spectrometer  due 
to  reflection  of  all  C2H4  gas  from  the  back  of  the  shutter.  At 
this  point,  the  shutter  is  removed  and  the  C2H4  signal  sharply 
decreases  due  to  adsorption  by  the  Si{100)  crystal.  The  region 
of  constant  PC9H4  due  to  an  adsorption  process  which  occurs 
with  a  constant  sticking  coefficient  equal  to  unity;  following 
this  region,  the  sticking  coefficient  drops  as  a  slower  adsorp¬ 
tion  process  occurs,  and  the  measured  PC2H4  signal  rises  asymp¬ 
totically  to  the  limit  shown  by  the  dashed  line.  At  this  point 
the  rate  of  adsorption  has  dropped  to  near  zero;  placing  the 
shutter  in  or  out  has  no  effect  on  the  measured  pressure  of  ethy¬ 
lene,  since  both  the  crystal  and  the  shutter  reflect  100%  of  the 
impinging  C2H4. 
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Kinetics  of  Adsorption  of  C2H4  by  Si(100) 


Schematic  Apparatus 


Beam  doser 

delivering  known  =CD 


c,h4  h„*=f 


Shielded  QMS 


Si(100)  crystal 
at  constant  T 


Shutter 


Schematic  Experiment 


Figure  2, 


Apparatus  for  Kinetic  Studies  of  Adsorption  on 
Si (100)-  Schematic  Experimental  Results. 


Given  that  the  sticking  coefficient  is  initially  unity,  it 
is  possible  to  determine  the  fraction  of  the  flux  of  C2H4  which 
intercepts  the  crystal,  using  the  fractional  initial  decrease  in 
PC2H4  uPon  opening  the  shutter.  In  our  case,  for  a  doser-to- 
cryscal  distance  of  1.0  cm,  the  fractional  interception  is  .61 
±  .03.  This  agrees  well  with  the  value  determined  from  calcula¬ 
tions  of  the  doser  performance  using  the  known  geometrical 
arrangement  (crystal  diameter  =1.6  cm;  capillary  array  doser 
diameter  =  1.06  cm),  where  a  fractional  interception  value  of 
0.52  ±  .05  is  calculated  [9], 


IV.  EXPERIMENTAL  RESULTS 

A.  Absolute  Coverage  of  C2H4  on  Si (100) 

A  number  of  measurements  have  been  made  of  the  saturation 
coverage  of  C2H4  on  Si (100).  These  are  listed  below  in  Table  I. 
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Table  I 


Saturation  Coverage  of  C2H4  on  Si (100) 
(Crystal  Area  =  1.96  cm2) 

Crystal  T  =  105  K. 


Experiment 

No. 

Adsorbate 

Number  of 
Molecules 
Incident 
(1014  cm”2) 

Number  of 
Molecules 
Adsorbed 
( 1014  cm”2) 

Saturation 
Coverage , 
C2H4/dimer 

1. 

472 

c2h4 

3.9 

1.21 

.36 

2. 

485 

c2h4 

2.5 

1.14 

.33 

3. 

491 

c2h4 

2.9 

1.27 

.  37 

4. 

492 

c2h4 

2.8 

1.20 

.35 

5. 

58. 

c2h4 

2.9 

1.29 

.38 

6. 

622 

c2h4 

2.4 

1.20 

.35 

7. 

623 

c2h4 

2.8 

1.21 

.36 

8. 

641 

c2h4 

2.8 

1.28 

.  38 

9. 

651 

c2h4 

2.2 

1.04 

.30 

10. 

653 

c2h4 

2.2 

1.04 

.31 

11. 

661 

c2h4 

2.3 

.94 

.28 

12. 

772 

c2d4 

2.  3 

1.09 

.32 

13. 

685 

c2d4 

4.1 

1.35 

.40 

14. 

687 

c2d4 

4.3 

1.07 

.31 

15. 

692 

C2d4 

4.3 

1.11 

.33 

16. 

695 

c2d4 

4.3 

1.10 

.32 

17. 

775 

c2d4 

4.0 

1.14 

.  33 

Average  =  0.34  ±  0.03 
C2H4/Si  dimer 
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B.  Adsorption  of  C2H4  via  a  Precursor  State  Mechanism 


The  initial  sticking  coefficient  of  unity  persists  until  a 
fraction  of  a _ monolayer  of  chemisorbed  ethylene  has  been 
adsorbed,  as  indicated  by  the  initial  flat  region  in  the  schema¬ 
tic  experimental  curve  shown  in  Figure  2.  This  may  be  understood 
from  Figure  3,  where,  in  the  right  hand  side  of  the  figure,  CoH^ 
migrates  over  filled  Si  sites  until  an  empty  chemisorption  site 
is  found  and  occupied.  Adsorption  via  a  mobile  precursor  state 
is  well  known  m  adsorption  studies  on  metals  and  has  been 
recently  reviewed  [10].  In  the  left  hand  side  of  Figure  3,  a 
schematic  potential  energy  diagram  is  shown  in  which  adsorption 
into  the  precursor  state  is  followed  by  two  possible  kinetic  pro- 
cesses:  ( 1 )_ desorption  from  the  precursor  with  rate  constant  kH ; 

and  (2)  chemisorption  onto  an  empty  site  with  rate  constant  kr. 


Detection  of  Mobile  Precursor  State 


Figure  3.  Mobile  Precursor  State  to  C2H4  Adsorption 
on  Si (100) . 
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Modeling  this  process,  we  obtain  the  relations: 


1  +  kr(l-er)  (1) 

where  Pac3  =  probability  of  adsorption 

F  =  flux  of  molecules  to  the  crystal 
kd»  kr  *  rate  constants  for  desorption  and  reaction 
8r  =  fraction  of  a  chemisorbed  monolayer 


^d  vd  ^d~^r  (2) 

kr  =  vr  e~  kT 

where  vr,  =  pre-exponential  factors,  and  Er,  E^  =  energy 
barrier  to  reaction  or  desorption. 

From  these  equations,  it  is  clear  that  we  can  find  the 
energy  difference  between  reaction  and  desorption,  (E^-Ej.)  and 
the  ratio  of  the  pre-exponential  factors  (v^/vf)  by  charac¬ 
terizing  the  different  probabilities  of  adsorption  with  varying 
crystal  temperature. 

The  rates  of  the  two  processes  will  vary  differently  with 
temperature  as  may  be  seen  from  the  different  activation  energies 
for  desorption  and  reaction  from  the  precursor,  and  this  fact  may 
be  used,  as  shown  in  equation  2,  to  determine  the  value  of  (E^  - 
Er).  In  our  experiments  we  have  used  a  temperature  range  from 
103  K  to  405  K  and  have  studied  the  influence  of  temperature  on 
the  initial  rate  of  adsorption  by  the  clean  crystal. 
Representative  data  are  given  in  the  top  of  figure  4  for  the 
various  temperatures;  an  Arrhenius  plot  yields  (E^  -  Er)  =2.9 
kcal/mole,  and  v^/vr  =27. 
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C2H4/Si(100)-Rate  of  Adsorption  versus  Temperature 


Arrhenius  Plot 


Figure  4.  C2H4  Adsorption  Kinetics  on  Si(100). 


If  we  accept  as  a  first  approximation  that  Ej  is  equal  to 
the  heat  of  vaporization  of  C2H4  =3.5  kcal/mole  [11],  then  Er  = 
0.6  kcal/mole.  To  the  extent  that  E^  is  greater  than  its  heat  of 
vaporization  because  of  additional  binding  in  the  second  layer, 
then  Er  will  also  increase. 
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Within  the  confines  of  the  transition  state  theory,  the 
ratio  of  preexponential  factors  has  physical  meaning.  Thus,  the 
entropy  of  activation  for  desorption  from  the  mobile  precursor  is 
greater  than  the  entropy  of  activation  for  reaction.  This 
implies  that  the  transition  complex  for  reaction  is  of  lower 
entropy  than  the  transition  complex  for  desorption,  that  is,  the 
reaction  with  the  adsorption  site  involves  a  transition  complex 
which  is  constrained  entropically  compared  to  the  transition 
complex  involved  in  desorption  from  the  mobile  precursor. 


C.  Desorption  of  C2H4  from  Si (100). 

Desorption  experiments  were  carried  out  using  both  ethylene 
and  deuterated  ethylene  (C2D4).  A  survey  of  all  desorbing  gases 
in  the  mass  range  2,  4,  12  -  64  amu  revealed  that  only  ethylene  and 
hydrogen  desorbed  from  the  surface  in  the  temperature  range  100  - 
1000  K.  Only  a  single  desorption  feature  is  observed  for  ethy¬ 
lene  as  shown  in  figure  5,  where  a  sharp  desorption  peak  with 
increasing  peak  temperature  in  the  range  554  K  to  591  K  is 
observed,  as  the  initial  coverage  is  increased.  The  desorption 
kinetics,  assuming  first  order  desorption,  was  fitted  according 
to  the  method  of  Chan  et  al  [12],  giving  (from  the  chemisorbed 
state)  =  38.0  kcal/mole  (figure  6),  and  =  10*3.7  sec”1,  in  the 
limit  of  zero  coverage. 
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THERMAL  DESORPTION  OF 
ETHYLENE  FROM  SILICON(IOO) 


300  400  500  600  700  800 

T(K) 


Figure  5.  Thermal  Desorption  of  C2H4  from  Si (100). 
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ETHYLENE  ON  SILICON(IOO) 


ENERGY  OF  DESORPTION  VS.  EXPOSURE 


Exposure  (10^  molecules/cm^) 


PRE-EXPONENTIAL  FACTOR  VS.  EXPOSURE 


Exposure  (10^  molecules/cm^) 

Figure  6.  Desorption  Rate  Parameters-  C2H4/Si ( 100 ) . 


D.  Decomposition  of  C2H4  on  Si (100). 

It  is  well  known  that  C2H4  will  inefficiently  decompose  on 
Si (100)  at  elevated  temperatures.  Thus,  in  the  work  of  Bozso  et 
al  [13],  an  efficiency  of  about  10“ 3  giC  units  per  C2H4  collision 
with  Si ( 100 )  at  940  K  was  measured.  Subsequent  work  by  Taylor, 
et  al  [141  remeasured  this  efficiency,  finding  an  efficiency  of 
about  10“’  -  10"5  SiC  units  per  C2H4  collision  at  970  K.  We  have 
shown  that  chemisorbed  C2D4  also  inefficiently  decomposes  upon 
heating  a  monolayer  above  the  desorption  temperature,  as  shown  in 
figure  7.  Here,  the  thermal  desorption  of  C2D4  was  interrupted 
at  various  temperatures  and  the  remaining  carbon  was  measured  by 
Auger  spectroscopy.  It  may  be  seen  that  near  600  K  the  carbon 
surface  concentration  drops  to  near  zero  as  C2D4  desorption  is 
completed.  Auxiliary  experiments,  in  which  the  amount  of 
desorbing  deuterium  is  measured  following  C2D4  desorption,  have 
shown  that  less  than  about  4%  of  the  C2D4  is  decomposed  on 
heating.  Such  low  levels  of  reactivity  could  well  be  attributed 
to  the  presence  of  surface  defects  on  the  Si (100)  crystal. 


Loss  of  Carbon  and  C2H4  Desorption 


versus  Temperature 


.w  a  = - W  u - u - w - W -  u. 

450  500  550  600  650  700  750 

TEMPERATURE(K) 


Figure  7.  Loss  of  Carbon  and  C2D4  Desorption  versus 
Temperature. 


V.  SUMMARY 

The  following  accomplishments  have  been  made  in  the  past 

year. 

1.  An  apparatus  for  quantitative  studies  of  the  kinetics  of 
adsorption  and  the  stoichiometry  of  adsorption  has  been 
constructed. 

2.  The  adsorption  of  ethylene  by  Si (100)  has  been  studied. 
It  has  been  found  that  the  adsorption  process  occurs  via 
the  involvement  of  a  mobile  precursor  state  which  ini¬ 
tially  causes  the  ethylene  sticking  coefficient  to  be 
unity  at  temperatures  near  100  K. 

3.  The  kinetics  of  desorption  and  reaction  from  the  mobile 
precursor  have  been  studied,  and  the  difference  in  acti¬ 
vation  energy  for  desorption  and  for  reaction,  (Ed  -  Er) 
has  been  measured  to  be  2.9  kcal/mole. 


4. 


The  kinetics  of  desorption  of  ethylene  from  its  chemi¬ 
sorbed  state  has  been  measured,  yielding  E^  =  38.0 
kcal/mole  and  =  10^.7  sec-*  for  first  order  kine¬ 


tics. 


5.  Nearly  100%  of  the  chemisorbed  ethylene  desorbs  from 
Si (100)  with  little  decomposition  to  surface  carbon 
occurring . 

6.  The  saturation  coverage  of  ethylene  is  estimated  to  be 
0.34  C2H4/Si  dimer.  Surface  models  consistent  with  this 
measurement  are  being  developed. 
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APPENDIX  IV 

CATHODOLUMINES CENCE  EXPERIMENTS 


I.  Description  of  the  Apparatus 


A  schematic  of  the  apparatus  is  shown  in  side  view  in  Figure  1. 
It  consists  of  an  electron  gun  for  excitation  of  the  emission,  UHV 
sample  chamber,  and  optics  to  spectroscopically  analyze  the 
cathodoluminescence  emission. 

Ultrahigh  vacuum  techniques  are  necessary  to  make  sure  that 

there  are  no  condensed  films  on  the  sample.  The  electron  gun  energies 

were  selected  to  excite  samples  as  thin  as  tens  of  angstroms,  and  almost 

all  gases  from  leaks,  desorption,  or  backstreaming  would  condense  on  the 

samples  that  are  maintained  at  cryogenic  temperatures.  Particularly, 

carbonaceous  pump  oils  irradiated  by  the  electron  beam  would  be  expected 

to  appear  very  "diamond-like.*  The  only  way  to  avoid  surface 

contamination  is  to  make  sure  that  monolayer  formation  times  are  longer 

than  experiment  times.  The  turbo-backed  ion  pump  and  sorption  pumps 

used  in  the  sample  chamber  maintain  pressures  of  5  x  10  ^  torr.  The 

cryogenic  sample  heat  shields  maintain  pressure  at  the  sample  at  10 
-12 

to  10  torr,  depending  on  temperature  of  operation.  The  turbo  is  also 
used  to  back  the  rotating  seal  on  the  sample  turret,  which  will  be 
discussed  later. 

The  monochromator  is  a  custom-made  1  meter  Ebert  vacuum 
instrument.  Fused  silica  and  sapphire  optics  permit  operation  at 
wavelengths  down  to  the  cutoff  in  air  at  about  1850  in  the  present 
configuration,  but  it  is  possible  to  vacuum-couple  the  instrument  to  the 
sample  chamber  if  required.  The  optical  arrangement  contains  a  beam 
splitter  which  superimposes  a  wavelength  calibration  line  spectrum, 
e.g.,  mercury,  on  the  measured  spectrum. 
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The  details  of  the  six-position,  rotatable  sample  turret  are 
shown  in  Figure  2.  One  of  the  six  sample  positions  is  occupied  by  a 
phosphor  which  is  used  for  alignment  by  positioning  a  laser  to  coincide 
with  the  electron  beam-excited  spot.  This  spot  is  focussed  on  the 
monochromator  entrance  slit,  so  when  samples  are  rotated  into  the  laser 
spot  they  are  automatically  in  the  correct  position  both  for  the 
electron  beam  and  for  the  analyzing  optics.  The  electron  beam  energy 
can  be  varied  from  300V  to  30  kV,  corresponding  to  a  wide  range  of 
stopping  distances  in  the  materials  of  interest  here. 

Three  samples  of  cathodoluminescence  spectra  measured  with  this 
apparatus  are  presented  below.  Figure  3  shows  the  spectrum  from  a 
natural  type  lib  diamond,  giving  a  broad  band  in  the  blue  with  few 
features.  The  sharp  lines  are  the  superimposed  calibration  lines,  and 
the  computerized  data  management  routine  interpolates  between  these 
lines  to  determine  the  wavelength  scale.  The  spectrum  shown  in  Figure  4 
was  taken  at  room  temperature  on  a  synthetic  polycrystalline  diamond 
film  deposited  at  Westinghouse  BAD  Center,  and  the  spectrum  in  Figure  5 
of  this  same  sample  taken  at  60K  shows  a  shift  of  the  main  band  toward 
higher  energies,  as  would  normally  be  expected. 

We  are  currently  analyzing  these  data  to  identify  the  bands  and 
features  of  the  spectra.  Also,  a  variety  of  samples  prepared  under 
tightly  characterized  conditions  is  being  obtained  in  order  to  study  the 
influence  of  dopants  and  defects,  and  to  provide  a  baseline  of  knowledge 
for  studying  the  materials  that  are  produced  on  this  program. 
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